Isotope effects in the Hubbard-Holstein model within dynamical mean-field theory 
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We study the isotope effects arising from the couphng of correlated electrons with dispersionless 
phonons by considering the Hubbard-Holstein model at half-filling within the dynamical mean- 
field theory. In particular we calculate the isotope effects on the quasi-particle spectral weight Z, 
the renormalized phonon frequency, and the static charge and spin susceptibilities. In the weakly 
correlated regime U /t < 1.5, where U is the Hubbard repulsion and t is the bare electron half- 
bandwidth, the physical properties are qualitatively similar to those characterizing the Holstein 
model in the absence of Coulomb repulsion, where the bipolaronic binding takes place at large 
electron-phonon couphng, and it reflects in divergent isotope responses. On the contrary in the 
strongly correlated regime U/t^ 1.5, where the bipolaronic metal-insulator transition becomes of 
first order, the isotope effects are bounded, suggesting that the first order transition is likely driven 
by an electronic mechanism, rather then by a lattice instability. These results point out how the 
isotope responses are extremely sensitive to phase boundaries and they may be used to characterize 
the competition between the electron-phonon couphng and the Hubbard repulsion. 



PACS numbers: 71.38.-k,63.20.Kr,71.38.Cn,71.38.Ht 

I. INTRODUCTION 

One of the most direct methods to establish the cou- 
phng of an electronic property with the dynamics of 
the underlying lattice is the isotope effect (IE), a a = 
—dhiA/dhiM, which probes the dependence of the elec- 
tronic quantity A on the ion mass M . For example, the 
observation of a large IE on the superconducting transi- 
tion temperature Tc has been an important evidence for 
a phonon-mediated pairing mechanism in conventional 
low-temperature superconductors. In the weak-coupling 
BCS model, the superconducting critical temperature is 
proportional to the Debye phonon frequency which scales 
as The resulting IE on ut, = 1/2, has been 

considered as one of the most important confirmation of 
the BCS theory 

The vanishingly small values of in optimally doped 
high-Tc copper oxides superconductors has initially in- 
duced a large part of the scientific community to be- 
lieve that the pairing mechanism in these compounds was 
mainly of electronic origin, and oriented the theoretical 
research towards purely electronic models, relegating the 
lattice to a secondary role. Later experimental investiga- 
tions showed however that small values of ar^ are pecu- 
liar of the optimal doping, whereas sizable values of ax^ , 
even larger than the BCS prediction, were reported in 
the underdoped phase. ^ Furthermore, the improvement 
of the experimental accuracy has allowed recently to es- 
tablish the presence of isotopic shifts also in several phys- 
ical quantities different from the superconducting Tc- 

Sizable isotope effects (lEs) have been indeed found 
on the zero temperature penetration depth Xi^ijd)^^ on 
the Knight shifty on the nuclear quadrupole resonanceji^ 



on the pseudogap temperature)^ and on the angular re- 
solved photoemission spectral These findings are par- 
ticularly remarkable since the Migdal-Eliashberg theory 
of the electron-phonon interaction predicts vanishingly 
small lEs other than ar^. For example, according to 
the BCS description of the superconducting state, the 
zero temperature penetration depth is given by Al(0) = 
(TO*/ns)^/2, where Ug is the density of the condensate 
and TO* is the effective electron mass. According to the 
Migdal-Eliashberg theory, to* = (1 -I- A)to, where to is 
the bare band mass and A is the electron-phonon cou- 
pling constant. Since A is independent of the ion mass, 
the isotope effect on the penetration depth is expected 
to be zero, in contrast with Refs. 2 3, 

The completely anomalous set of isotope dependences 
signals a relevant electron-phonon coupling, but at the 
same time highlights that the interplay with strong cor- 
relations and possibly with other features of the cuprates 
needs to be taken into account. Different theoretical 
models have been proposed to account for these anoma- 
lous lEsi^i^i^^i^-^i^^i^'^i^^ The analysis has been mainly 
concentrated on purely electron-phonon narrow-band 
systems, as the Holstein model which however neglects 
possible effect of the anisotropy of the electron-phonon 
interaction considered in Ref. 15. 

The lEs on to,* in interacting electron-phonon systems 
in the absence of electron-electron interactions have been 
investigated in details in Refs. 8.9- 10. More recently, the 
ability of the lEs in revealing different physical regimes 
of electron-phonon systems has been discussed in Ref. 

where remarkable anomalous lEs on to* and on the 
dressed phonon frequency Q.q were reported as the sys- 
tem enters in a polaronic regime. In order to focus on 
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the metallic properties and to clarify the origin of these 
anomalous lEs, however, a spinless Holstein model in the 
DMFT approximation was considered there, which en- 
forces the metallic character in the whole space of pa- 
rameters. 

In the opposite case of strong correlations the electron- 
phonon interaction in the weak-intermediate limit has 
been treated for instance by means of an adiabatic 
theory on small clustersji^ or of coherent potential 
approximation.^^ In this case the IE reflects in a smooth 
modification of high energy features such as the Hub- 
bard bands. In the strongly correlated antiferromagnetic 
phase lEs on the single hole spectral properties have been 
also analyzed with diagrammatic quantum Monte Carlo 
techniques 

The scope of this paper is to provide a nonperturba- 
tive analysis of the isotope effects on different observ- 
ables in a regime in which strong correlation and inter- 
mediate/strong electron- phonon interaction coexist. To 
this aim we employ the dynamical mean-field technique 
(DMFT), which is powerful tool to investigate the non- 
perturbative regimes of strongly correlated and electron- 
phonon systems. In particular, in this paper we consider 
the metallic regime of the Hubbard-Holstein model, de- 
liberately excluding broken-symmetry phases. We also 
limit ourselves to the half-filled systems, where in DMFT 
the correlation-driven Mott transition or the electron- 
phonon driven pairing transition take place as one or the 
other coupling becomes large, and we follow the behav- 
ior of the isotope effects when these phase boundaries 
are approached. It is well known that DMFT becomes 
exact only in the limit of infinite coordination number, 
and it has to be viewed as an approximate method in 
finite dimensions. Yet, DMFT allows us to access the 
full local quantum dynamics, and to deal with many dif- 
ferent energy scales. Moreover the assumption of a local 
self-energy, implicit in the DMFT, makes this technique 
not suitable to investigate the momentum dependence 
of the physical quantities. This is for instance the case 
of Frolich-like electron-phonon Hamiltonians and of the 
long-range Coulomb interaction, which gives rise to im- 
portant screening effects here neglectedi^ This is also 
the case of the momentum dependence of the electron- 
phonon interaction induced by the electronic correlation 
j^ggjgij22j23^ which cannot be properly take into account 
within the context of the DMFT. 

A general belief is that the electronic correlation com- 
petes with the electron-phonon interaction leading to a 
reduction of the phonon effects on the electronic proper- 
ties with a consequent decrease of the magnitude of the 
possible lEs. We show that in particular regimes, close 
to phase boundaries, the Hubbard repulsion can actually 
enhance the dependence of many physical quantities on 
the phonon frequency, yielding an increase of the lEs. 
Indeed the isotope effects turns out to be extremely sen- 
sitive to the closeness of metal-insulator transitions, and 
emerge as a tool to investigate the physics underlying the 
different regimes of the phase diagram. 



This paper is organized as follows. In Sec. |n]we in- 
troduce the Hubbard-Holstein model and the dynamical 
mean field method. In Sec. IIIII and Sec. |^ we report 
our results for the weakly correlated and strongly corre- 
lated regimes, respectively, while a study of the lEs as 
a function of the Hubbard-[/ repulsion is the subject of 
Sec. The conclusions are summarized in Sec. IVII 



II. OUR APPROACH 

The Hubbard-Holstein model is defined by the Hamil- 
tonian: 

+9 X! "io- i^^i + "l) + ^0 X! 
ia i 

where t is the hopping amplitude between nearest neigh- 
boring sites, U is the Hubbard electron-electron repul- 
sion, g is the electron-phonon matrix element, and ojq 
the Einstein phonon frequency. 

We can identify three independent dimcnsionlcss pa- 
rameters, the electron-phonon coupling constant A = 
2(7^/(wot), the adiabatic ratio 7 = oj^/t, which 
parametrizes the role of the lattice quantum fluctuations, 
and the dimensionless Hubbard repulsion U /t. It is im- 
portant to note that ujq oc I/^/M and g cx l/v^McJo, 
where M is the ionic mass, so that the only physical 
quantity which depends on M is the adiabatic ratio 7. 
Finite lEs stem thus from an explicit dependence on 7 
and they are hence a direct signal of the role of the lat- 
tice quantum fluctuations. The locality of the interac- 
tion terms makes the Hubbard-Holstein model particu- 
larly suitable to be investigated by dynamical mean-field 
theory.^'* In this approach the spatial fluctuations are 
frozen, i.e., all the lattice sites are assumed to be equiv- 
alent, but the local electron and phonon dynamics are 
fully accounted for. In this sense the method generalizes 
to a quantum framework the classical mean-field theories, 
and it becomes exact in the infinite coordination limit. 
Hence the lattice model is mapped onto a quantum im- 
purity model in which an interacting site is embedded 
into a non-interacting bath, whose spectral function has 
to be self-consistently determined. In our case, the impu- 
rity site displays both Coulomb repulsion and a coupling 
to a local phonon. The self-consistency equation con- 
tains all the information about the original lattice model 
through the bare density of states. As often, here we 
consider an infinite coordination Bethe lattice with half- 
bandwidth t,^^ which reproduces the finite bandwidth of 
finite-dimensional systems. For the normal paramagnetic 
phase, the self-consistency condition on the Bethe lattice 
reads 

AH = ^G(^), (2) 
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where G is the Green's function on the impurity site, 
which corresponds to the local component of the lat- 
tice Green's function, and A is the bath hybridization 
function. Within DMFT the electron and phonon self- 
energies are local, i.e., momentum independent. 

To solve the impurity model we use the exact diagonal- 
ization methodj-fi in which the dynamical DMFT bath is 
described in terms of a finite set Ns of Anderson impu- 
rity levels, while a finite cut-off is imposed on the total 
number Npi^ of the phononic Hilbert space. Typical val- 
ues we considered are Ng = 9, 10, and A'ph ~ 20. Both 
these quantities were varied to check the effective sta- 
bility of the results. The self-consistent DMFT set of 
equations was solved in the Matsubara space where a fic- 
titious temperature T plays the role of an energy cut-off. 
Particularly small value of T up to T/t ^ 1/400 were 
needed to ensure a high accuracy of our results. 

In this paper we investigate unconventional lEs on 
both single-particle and two-particle quantities. In par- 
ticular, concerning the one-particle properties, we focus 
on the quasi-particle spectral weight Z, which in DMFT 
coincides with the inverse of the effective electronic mass 
m*/m due to the momentum independence of the self- 
energy E 

and on the renormalizcd phonon frequency Qq obtained 
as21 

/r!oy_ ^0) 

V Wo / t^o 

In the above expressions, E(w„) is the electronic self- 
energy and D(ujm) is the phonon propagator D{uj„i) — 
— {x{u!m)x{—uJm)), whcrc X = (a -I- a^) is the dimcn- 
sionless lattice coordinate operator, ujn and a;,„ are 
respectively fermionic w„ = 7rT(2n -I- 1) and bosonic 
LUjn — 27rTm Matsubara frequencies. In practice Z is 
computed by linearly extrapolating from the first Mat- 
subara frequency. 

The knowledge of the renormalizcd phonon frequency 
51o allows us to evaluate also the static component of the 
local charge susceptibility p = p{uJm = 0), where p{ujm) = 
{n{LiJm)n{—uJra)) , n being the density operator on the im- 
purity. It is indeed easy to prove that /9 = [l—ujU^D/Xt. 
We evaluate also the local static spin susceptibility x = 
x(t^m ~ 0) for the spin response function x(^ni) = 

{S:^{uJra)Sz{-i^rn)), whcrC S^, = (1/2) X]a/3 C^CT^^C^ iS 

the z-component of the local electron spin (a^ being the 
Pauli matrix fa). Finally, we consider also the local lat- 
tice probability distribution function P{x) = |(0o|a;)p, 
where \x) is the eigenstate of the lattice coordinate op- 
erator X. In our exact diagonalization approach we can 
evaluate P{x) = J2n,m'^nix)ipTnix){4>o\n){m\(l)o), where 
|0o) is the ground state, \n) are the eigenstates of the 
harmonic oscillator and ipn{x) the corresponding eigen- 
f unctions. 



Despite its simplicity, the Hubbard-Holstein model has 
quite a rich phase diagram, in which different phases and 
regimes appear. In order to have a better comparison 
between the different regimes, we organize our discussion 
according to degree of correlation. In the next section 
we focus on the weakly correlated regime of the Holstein- 
Hubbard model, where the presence of the on-site Hub- 
bard repulsion changes only quantitatively the results of 
the pure Holstein spinful model, but it does not give rise 
to new features. We then discuss in section llVl the regime 
of strong electronic correlation where the interplay be- 
tween electron-phonon and electron-electron interaction 
induces qualitative changes on the behavior of the lEs. 



III. WEAKLY CORRELATED REGIME 

Already in the absence of correlation, the introduc- 
tion of the spin degeneracy introduces additional physics 
with respect to the spinless system, whose isotope effects 
have been studied in Ref. [l^ In the spinless case in- 
deed the polaron crossover is not associated with a metal- 
insulator transition (MIT) which is only asymptotically 
approached for A — s- cx). The persistence of a metallic 
character in the polaronic phase is reflected in the finite 
value of both the quasi-particle weight and the renormal- 
izcd phonon frequency.^^ On the other hand in the spin- 
ful case the residual attractive interaction between the 
electronic or polaronic charges leads to the possibility to 
form a insulating phase of localized pairs which in the 
adiabatic regime (7 <C 1) are accompanied by strong lat- 
tice deformation (bipolarons)^ This occurs within the 
DMFT approximation for A larger than a finite critical 
value where the quasi-particle weight vanishes at Amit 
signaling a pair-MIT j^Si In the adiabatic limit (7 = 0) 
the pair-MIT has a precursor in the phonon softening 
which occurs at Aph-'^" As the bare phonon frequency 
is increased, the phonon softening tends to occur more 
closely to the pair-MIT. 

In Fig. ^we show the quasi-particle spectral weight Z, 
the renormalizcd phonon frequency CIq, the local static 
spin susceptibility x ^J^d the local static charge suscep- 
tibility p as functions of the electron-phonon coupling A 
for different weak-intermediate values of the Hubbard re- 
pulsion U and for adiabatic ratios 7 = 0.1, 1.0. xo and po 
represent respectively the local spin and charge suscepti- 
bility in the absence of electron-phonon and Hubbard in- 
teraction. As mentioned above, the main difference here 
with respect to the spinless Holstein model is the pos- 
sibility to form local bipolaronic pairs of electrons with 
opposite spin. This leads to a metal-insulator transition 
at a finite electron-phonon coupling Amit where Z — > 0. 
The metal-insulator transition is accompanied by a soft- 
ening of the phonon frequency fio — > for a value Aph 
of the electron-phonon very close to Amit , followed by a 
sudden increase for A > Xp^m^iSL^ 

The approach to the bipolaronic metal-insulator tran- 
sition is reflected also in the two-particle quantities x 
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FIG. 1: (color online) Behavior of the quasi-particle spectral 
weight Z, the renormalized phonon frequency Slo, the local 
static charge susceptibility p and the local static spin suscep- 
tibility X s-s functions of the electron-phonon coupling A for 
different weak-intermediate values of the Ifubbard repulsion 
U . In the left panel the curves are plotted for the adiabatic 
ratio 7 = 0.1 while in the right panel the adiabatic ratio is 
7 = 1.0. 

and p. In particular, the increase in the number of sin- 
glet local electron pairs close to the transition leads to 
a suppression of the local spin susceptibility and to an 
enhancement of the local charge susceptibility, so that x 
has a sharp drop whereas p nearly diverges at Aph. 

In this context the Hubbard repulsion gives rise to a 
further reduction of metallic properties at small A (de- 
crease of Z) while, on the other hand, it disfavors double 
occupancy leading to a shift to larger A of the critical 
values of the el-ph coupling and to an increase (decrease) 
of the spin (charge) susceptibility. Finally, increasing the 
quantum lattice fluctuations (7 = 1.0, right panels) leads 
to a further increase of the critical electron-phonon cou- 
plings. 

In Fig. 121 we show the corresponding isotope coeffi- 
cients 

ax = -AlogX/AlogM = (l/2)AlogX/Alog7, (5) 

where X is Z, fioi X ^'^d p. We have chosen to compute 
the lEs using a finite value /S.^/^ = 0.15. 

The behavior of lEs approaching Amit and Aph is 
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FIG. 2: (color online) Plot of the isotope coefficients for the 
same quantities defined in Fig. 



depicted in Fig. |21 The most striking effect is the 
near divergence of the isotope coefficients as the sys- 
tem approaches the pairing metal-insulator transition. 
As a matter of fact, while az has a true divergence as 
A Amit J for the isotope effects on the other quantities 
the divergence is cut-off by the not complete suppression 
of charge fluctuations. From the practical point of view 
however both the MIT and bipolaronic instabilities are 
accompanied by a huge increase of the magnitude of the 
anomalous isotope effects. 

The key point to understand this issue is that an iso- 
topic shift of is reflected in a finite shift to higher A of 
both the bipolaronic MIT (Amit) and the phonon soft- 
ening (Aph). It obviously descends from Eq. ^ that the 
divergence of a quantity is accompanied by a divergence 
to —00 of the corresponding isotope coefficient a, while 
a positively divergent isotope coefficient is associated to 
a vanishing observable. According this view we can un- 
derstand the behaviors of Fig. [3 Namely we find that 
az ^ 00 as Z ^ at Amit and a-^ nearly diverges as x 
rapidly drops at Amit- In similar way the near vanishing 
of r^o and the near divergence of p, at Aph are reflected 
in Q!o„ — !■ -|-oo and ap — > —00. The following drop of a^^ 
and the jump of ap for A > Aph reflect the sharp behavior 
shown in Fig^of the corresponding quantities. 

Increasing U leads to an increase of Amit, while in- 



5 




"0 0.1 0.2 0.3 0.4 0.5 

X 

FIG. 3: (color online) Behavior of single occupation number 
P ~ \{t I0o)|^ + |(i \4>o)\'^ as function of A in the weak el-ph 
coupling regime (A ~ 0.5) for different adiabatic ratio 7 and 
for [7 = 0. 



creasing luq determines an overall smoothing of these fea- 
tures. An interesting point to be noted here is the non 
monotonic behavior of the IE on the spin susceptibility x, 
which shows a negative sign at small A before approach- 
ing a near divergence towards -f 00 close to the metal- 
insulator transition. On the contrary ap remains nega- 
tive up to Aph. The negative sign of and ap for small 
A can be understood in a perturbative framework. The 
electron-phonon interaction introduces a positive Stoner- 
like contribution to the effective spin susceptibility: 



X 



Xo 



1 - /el-el + l!" 



(6) 



cl — ph 



where /oi-ci oc U arises from the electron-electron repul- 
sion and the electron-phonon term I^^_pYi 7A takes into 
account the explicitly isotope dependence due to the fi- 
nite bandwidth and vertex corrections arising from the 
electron-phonon interactioniSiS In the same RPA fashion 
the charge susceptibility reads 



P = 



Po 



l + L 



el— el 



^ol-ph 



(7) 



where I^i_pi^ — —^ + Iei-ph- Thus at a perturbative level 
the only source of IE in p and x comes from I^i_pi^, which 
can be shown to be the same as in both quantities leading 
to the same negative IE. 

The sign of the isotope coefficient on the spin suscep- 
tibility reflects thus two different physical regimes where 
lattice fluctuations triggered by ojq play a different role. 
Namely, in the perturbative small A regime lattice fluc- 
tuations disfavor single spin occupancy (Fig. |3J), leading 
to a decrease of x a-nd to a negative isotope coefficient. 
On the other hand the physics in the A ^ Amit region 
is dominated by the proximity to the bipolaronic MIT. 
In this context lattice fluctuations relax the phonon- 
induced electron trapping favoring the unpaired spin sin- 
glet states. Along this line the sign of the isotope coef- 
ficient on X can be thus used as a useful trademark to 
characterized the different physical regimes. The ability 
of the IE to characterize the physical properties of differ- 



ent regimes will be even more remarkable in the strongly 
correlated regime. 



IV. STRONGLY CORRELATED REGIME 

In the previous section we have seen that for small U 
the physical properties of the system are qualitatively 
similar to those of the pure Holstein model, where turn- 
ing on a small Coulomb repulsion leads essentially only to 
an increase of the critical value of A for bipolaron pairing. 
The situation is significantly changed in the strongly cor- 
related regime in which the Hubbard repulsion exceeds 
the electron-phonon coupling. In this case a more rea- 
sonable starting point is the physics of the pure Hubbard 
model whose metallic region, close to the Mott transi- 
tion, can be described in terms of a narrow resonance 
at Fermi level accompanied by the two high-energy Hub- 
bard subbands. In this framework the electron-phonon 
interaction has two opposite effects: on one hand it pro- 
vides a further source of scattering which would decrease 
the quasi-particle scattering weight Z; on the other hand 
it acts as a effective attraction between the electrons, 
competing with the Hubbard repulsion;-- In the strongly 
correlated regime, when the reduction of the effective re- 
pulsion prevails over the electron-phonon renormalization 
of the quasi-particle properties, increasing A leads to an 
enhancement of Zm^i^ 

Another interesting feature of the Hubbard-Holstein 
model pointed out in Ref. |^ in the strongly correlated 
regime is the change of the order of the bipolaronic MIT 
at 7 = 0.1. The bipolaronic transition, which is contin- 
uous (second order) for U/t ^ 1.5, becomes indeed dis- 
continuous (first order) for L/ 1.5. In this latter regime 
all the observable have thus a jump at Amit- The evo- 
lution from a second order to a first order transition has 
been discussed in terms of a Landau-Ginzburg theory 
where the first/second order transition is driven by the 
local lattice configuration. From a numerical point of 
view, however, it is hard to distinguish, on the basis of 
observables like Z or ^Iq a first order transition from a 
sharp second order one. We show that the analysis of 
the phonon IE is a powerful tool to distinguish in a very 
clear way the first order or the second order character of 
the bipolaronic MIT. This is even more interesting since, 
as we are going to show, the first order transition is not 
related to a phonon mechanism while it is probably of a 
electronic origin. 

In Fig. ^ we show Z, Qq, x Si,iid p as function of A 
in the regime of strong correlation for 7 = 0.1 and for 
7 = 1.0. As mentioned above, in this strongly corre- 
lated regime, increasing A, as long as the system does 
not reach a bipolaronic regime, results in a reduction of 
the effective Hubbard repulsion and in an increase of Z. 
It is however interesting to note that the behavior of the 
other quantities, namely JIqi X a-nd p, as function of A 
and U does not change qualitatively with respect to the 
weakly correlated case. Concerning the character of the 
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FIG. 4: (color online) Behavior of the quasi-particle spectral 
weight Z, of the renormalized phonon frequency flo and of the 
charge and spin susceptibilities, p and as functions of the 
electron-phonon coupling A for different intermediate-strong 
values of the Hubbard repulsion U. In the left panel the curves 
are plotted for the adiabatic ratio 7 = 0.1 while in the right 
panel the adiabatic ratio is 7 = 1.0. 



bipolaronic MIT, we find that the transition is always of 
the second order for any U a.t j — 1.0. On the other 
hand, due to the lack of precursor effects in the metallic 
phase, it is hard to predict the character of the transition 
at 7 = 0.1. 

The corresponding lEs on the quantities Z, flo, x ^-nd 
p is reported in Fig. [S| We note a qualitatively different 
behavior of the isotope coefficients in the strongly cor- 
related regime at 7 — 0.1 as compared with the weak 
U regime and (partiaUy) with the case 7 = 1.0. Let 
us discuss for the moment the highly correlated cases 
U/t — 2.0, 2.5. Most striking is the behavior of the 
phonon and of the charge susceptibility isotope coeffi- 
cient which shows a trend diametrically opposite to the 
weak-C/ case, namely a marked reduction (increase) of 
Q^sio (ctp) to be compared with the divergences a^o — >■ 00, 
Q!p — > ^00 for small U. Among the other anomalous fea- 
tures we note a typical downturn of az as function of A, 
which is less marked as U increases. On the contrary, 
the upturn of the isotope coefficient on the spin suscep- 
tibility which is evident at small U is here absent. Quite 
interesting is also the U/t = 1.5 case which is very close 




FIG. 5: (color online) Plot of the isotope coefficients for the 
same quantities reported in Fig. 0] 



to the value where the bipolaron MIT changes from the 
second order to the first order. In this case the lEs show 
similar trends {aao < 0.5, ap > 0, downturn of az) as 
for the first order cases in a large range of A, with a fi- 
nal upturn (downturn for ap) towards the characteristic 
divergences of the second order. This trend is clearly vis- 
ible also at 7 = 1.0 (although the bipolaronic transition 
is always of the second order) pointing out an partial un- 
derlying tendency towards a first order instability. The 
larger amount of the lattice fiuctuations in this case as 
compared with 7 = 0.1 prevents however the occurrence 
of the first order instability and enforces the second order 
character. 

As we are going to show, the physical origin of these 
anomalous behaviors can be related to different phonon 
regimes, although they do not rule directly the change 
between the first and second order character of the tran- 
sition. In order to clarify in more details this issue we 
evaluate the local lattice probability distribution function 
(PDF) P{x) which provides information about the lattice 
wavefunction of the ground state. In Fig. El (upper panel) 
we show the PDF for two representative cases charac- 
terized by significant lattice distortions induced by the 
electron-phonon interaction, resulting in a non-gaussian 
P{x). In the strictly adiabatic case, polaronic effects for 
instance are refiected by a bimodal shape of P{x), where 
the two maxima are associated to different electron oc- 
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FIG. 6: (color online) PDF lattice function P(x) and lattice 
potential V(x). Solid black line: U jt = 2.0, A = 2.0, dashed 
red line: U jt = 1.5, A = 1.7. 



cupancies (e.g., in spinless systems n = and n = 1, 
in spinful systems n = and n — T). The broadness of 
these structures induced by quantum lattice fluctuations 
however can make hard to resolve the effective extent of 
the polaronic lattice distortions. To overcome this prob- 
lem we extract an "effective" lattice potential V(x) by 
solving the inverse Schrodinger problem 



1 

^ 

' 2M 



V{x) 



(8) 



where '0(x) — P{x) and M is the atomic mass. Fig. 
shows the lattice distribution function P[x) and the cor- 
responding lattice potential Vix) for two representative 
cases where the multi- valley structure of Vix) is not re- 
flected in the shape of P{x). We stress that the above de- 
fined Vix) reproduces the physical lattice potential only 
in the strict adiabatic limit 7 — > 0. It can be consid- 
ered however to be still a valid approximation at finite 
and small 7 providing thus useful insights for the true 
(unknown) lattice potential. 

In Fig. [7|we report the evolution of the lattice potential 
V{x) as function of A for three representative cases of the 
Hubbard repulsion, U jt — 1.0, 1.5,2.0, and for 7 = 0.1. 
All the shown A's belong to the metallic phase. The 
largest value of A, in each panel of the left side of Fig. [7| 
corresponds to the MIT. 

The first case (panel a) shows, in agreement with a 
second-order transition scenario, a continuous flattening 
of the lattice potential as A increases. This picture shows 
in the clearest way the lattice instability at Aph ~ 1.3, 
where the quadratic term of V(x) ~ aix^ -t- 042;* -I- . . . 
vanishes and the renormalized phonon frequency f^o 0- 
For Aph < A < Amit the systems is still metallic and 
a double-well structure of the lattice potential is devel- 
oping. This corresponds to the phonon hardening back 
shown in Fig. ^ The electronic metal-insulator tran- 
sition occurs for Amit > Aph when the double-well po- 
tential is already established and the MIT is caused by 



> 




FIG. 7: (color online) Lattice potential V(x) in differ- 
ent physical regimes. Upper panels: U jt = 1.0, A = 
1.0, 1.1, 1.2, 1.3, 1.34; (right: A = 1.34). Middle panels: U jt = 
1.5, A = 1.4,1.5,1.6,1.7,1.71246; (right: A = 1.4, , 1.71246). 
Lower panels: U jt = 2.0, A = 1.6, 1.7, 1.8, 1.9, 2.0; (right: 
A = 2, 0). The horizontal lines in panels b,d,f marks the zero 
point lattice energy for the corresponding lattice potential. 



the vanishing of coherent tunneling between the two min- 
ima. A crucial role is played by the zero point phonon 
energy wo/2 (horizontal lines in panel b) which favors 
tunneling processes. In particular, in the latter situation 
where the double- well is already formed, an increase of 
tjQ — ti^o + ^1^0 J due to an isotope shift, further enhances 
the metallic character, enhancing Z and resulting in a 
positive isotope coefficient on Z . At the same time the 
system, increasing wo, will probe a less shallow potential 
(panel b) resulting in a stronger hardening of SIq than 
the harmonic one and leading to > 0.5. 

A more complex situation is found for U jt = 1.5. Here, 
more than a softening of the quadratic term 022;^, in- 
creasing A the lattice potential V{x) presents two flat 
regions at finite x which evolve first in two relative min- 
ima and then in two absolute minima lower than the 
central one (panel c). In this framework we distinguish 
two regimes. The first one is when the zero point en- 
ergy is roughly lower than the flat regions (or minima) 
of the potential (A < 1.6); in this situation increasing 
luq makes the lattice wave-function more distorted, the 
metallic character is reduced and the phonon energy is 
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softened with respect to the harmonic one. This phe- 
nomenology is thus reflected in the downturn of az and 
of < 0.5. A second regime is reached for A closer 
to the MIT (A ^ 1.6): here the minima are lower than 
the zero point energy, a phenomenology similar to the 
double- well potential is restored and az > 0, a^o > 0.5. 

The results for U jt — 2.0 (panels e,f) are essentially 
similar to those for U jt — 1.5 for small A. Here however 
the first-order metal-insulator transition occurs well be- 
fore the relative minima at finite x reach the range of the 
zero point energy. As a consequence the previous con- 
siderations apply, namely az shows an downturn (but it 
still remains finite) and < 0.5. 

The evolution of the lattice potential as function of U 
and A sheds interesting light on the origin of the first- 
order/second-order character of the MIT. It has been 
indeed claimed that both the second-order at small U 
and the first-order transition at large U could be related 
to a phonon-driven mechanism^ In a phonon-driven 
Ginzburg-Landau theory the character of the transition is 
thus determined by the dependence of the free energy on 
the lattice coordinate, F{x) ~ a2x'^ /2 + a4X* /A + qqx^ /Q, 
which in the adiabatic limit reduces to the lattice poten- 
tial V{x). The second-order transition at small U is thus 
associated with the development of a double-well poten- 
tial which occurs for a2 < and ai,ae > 0, whereas a 
first order one is expected to correspond to a three-well 
structure with absolute minima at xq ^ 0, condition re- 
alized for 02, ag > 0, 04 < —^/Aa^OQ. Our results shown 
in Fig. 13 confute however this picture. 

At U ^ 1.5 the second order MIT transition occurs af- 
ter the development of the double well in V{x). MIT 
here is indeed a by-product of the lattice distortions as 
was shown by adiabatic analysis of Ref. and is due 
to the vanishing of the coherent tunneling processes be- 
tween the minima of two double-well lattice potential. 

For [/ 1.5 instead the developing of the three- well 
shape of the lattice potential is not sufficient to explain 
the first-order character since the MIT occurs far be- 
fore the local minima with finite lattice distortions at 
X become absolute minima. Moreover the energy 
of the local minima at the transition become higher and 
higher in energy as U increases, precisely where the first- 
order character becomes more marked. These considera- 
tions suggest thus that the origin of the first-order metal- 
insulator transition in the strongly correlated regime is 
likely ascribed to an electronic origin more than some 
phonon-driven instability. An intrinsic difference be- 
tween the weakly and the strongly correlated metallic 
phase is the presence of a spin-singlet character in the 
metallic case pointed out by the presence of the Hubbard 
bands coexisting with the Kondo peak in the strongly 
correlated regime. If the metallic state has a pronounced 
spin-singlet character the transition to a bipolaronic pair 
state (which is instead a charge pseudo-spin singlet) can 
be viewed as a level crossing between this two different 
electronic states therefore leading to a first order transi- 
tion. 





FIG. 8: (color online) Left panel: quasi-particle spectral 
weight Z and renormalized phonon frequency Qo as functions 
of Hubbard repulsion U for two different A — 0.7 (empty 
diamonds) and A = 1.5 (empty circles). Right panel: the 
corresponding lEs az and aoo- AH curves are plotted for 
ujo/t = 0.1. In the right-top panel are also shown the ana- 
lytical behavior predicted by Eq. Q (dashed lines) by using, 
respectively, Uc2/t = 2.87 for A = 0.7 and Uc2/t = 2.90 for 
A = 1.5. 



V. [/-DEPENDENCE 

In the previous section we have analyzed in details the 
behavior of the lEs on the physical quantities as function 
of the electron-phonon coupling A. We have shown in 
particular that the lEs are quite sensitive to the different 
regimes of the electronic correlations, and we have dis- 
tinguished the weak from the highly correlated regime. 
The aim of this section is to address this issue in more 
details by studying the dependence of the lEs on the 
Hubbard repulsion U. We shall focus on the IE on the 
quasi-particle spectral weight Z and on the renormalized 
phonon frequency f2o which show the most characteristic 
features. We consider two values of A, representative of 
two different physical regimes; A = 0.7, for which the 
system is far from the bipolaronic phase even in the ab- 
sence of any Hubbard repulsion; and A — 1.5 for which 
value instead the ground state is a bipolaronic insula- 
tor at J7 = 0. In this case one needs a finite value of 
the Hubbard repulsion {U ^ 1.25) to contrast the el-ph 
driven attraction and to restore the metallic character. 
Both the cases of course undergo a Mott metal-insulator 
transition as U approaches some critical value Uc- 

In Fig. |Hlwe report the evolution of the quasi-particle 
spectral weight Z and of the renormalized phonon fre- 
quency fio as functions of U for the two chosen A. Let us 
discuss first the behavior of flo, which shows a monotonic 
behavior for A = 0.7. This trend is representative for the 
generic case of weak electron-phonon coupling: due to 
the suppression of the double/empty sites, increasing U 
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leads to a reduction of the charge fluctuations and to 
the consequent enhancement of the renormahzed phonon 
frequency which results less screened. More intrigu- 
ing is the analysis of Z which, unlike the pure Hubbard 
model, shows a non-monotonic behavior as function of 
IJ . In the previous section we have discussed how, in 
the strongly correlated regime close to the Mott MIT, 
switching on the electron-phonon interaction leads to an 
increase of Z. This trend has been explained as a result 
of the competition between the attractive el-ph coupling 
and the repulsive Hubbard termi2ii24 Fig. |H1 shows that 
a similar mechanism applies also close to the bipolaron 
MIT where the phonon-driven attraction is the strongest 
interaction, and the Hubbard repulsion partially counter- 
acts it, so that Z increases with the repulsion for small 
\J . For larger U however the repulsion overcomes the at- 
traction, the system approaches the Mott MIT and Z re- 
covers the standard behavior (Z ^ 0). These trends are 
also more marked for A = 1.5 where, as mentioned above, 
the system is a (bipolaronic) insulator in the absence of 
Hubbard repulsion, and a finite f/ ^ 1.2 is required to 
have a metallic state. Most evident is of course the com- 
petition between U and the electron-phonon interaction, 
which becomes a sharp transition from Z = for J7 ^ 1.2 
to a maximum value Z ~ 0.5 for U ~ 1.5. 

In the right panels of Fig. |Slwe show the corresponding 
lEs az, afio- Also here the behavior of the isotope effect 
on the quasi-particle spectral weight is particularly inter- 
esting. In section IIIII we have shown that giant lEs on 
the quasi-particle spectral weight Z are expected as the 
system approaches the bipolaronic metal-insulator tran- 
sition. From an intuitive point of view we can expect 
that the electronic correlation, driving the systems far 
from the bipolaronic transition, reduces the magnitude 
of the anomalous IE az- This picture is confirmed in 
Fig. IHl where we see that az is strongly reduced as soon 
as U moves the system away from the bipolaronic MIT. 
This trend is however quite soon followed by a further 
enhancement of the anomalous IE on Z which becomes 
larger as the systems approaches the Mott transition, re- 
vealing that the strongly correlated regime close to the 
Mott metal-insulator transition is also highly sensitive to 
the phonon dynamics. This counterintuitive result can 
be qualitatively understood by following the analysis of 
Ref . 13 which showed that the strongly correlated regime 
close to the Mott MIT in the 'presence of electron-phonon 
interaction could be described in terms of a pure Hub- 
bard model with a properly rescaled Hubbard repulsion 
C/eff = U — rjXt, where ij = 2ujo/{2ujo + U). Since the 
quasi-particle spectral weight Z is known to scale lin- 
early with Uc2 - C/ in DMFT for U ^ U2c: where Uc2 
is the critical value of the Hubbard repulsion at which a 
Mott metal- insulator transition occurs, we obtain thus 

az = (9) 

This expression shows that the isotope effect on Z di- 
verges as l/Z as UcS Uc2^^^ The analytical behavior 



predicted by Eq. |^ is also shown in the right-top panel 
of Fig. |S1 (dashed lines). The critical values Uc2 for each 
A have been chosen in order to have the best fit with 
the DMFT numerical results: namely Uc2/t = 2.87 for 
A = 0.7 and Uc2/t = 2.90 for A = 1.5. 

The crossover between the two different correlation 
regimes is also pointed out by the analysis of the behav- 
ior of the lEs on the renormalized phonon frequency ao^ . 
Most evident is the case of A — 1.5 where in the absence 
of Hubbard repulsion the ground state is a bipolaronic in- 
sulator. Here the non-monotonic behavior of a^a reflects 
two different regimes: i) a flrst sharp region U ^ 1.36 
where a^a ^ 0.5. In this regime the system is character- 
ized by a double-well lattice potential with a zero point 
energy larger than the barrier between the two minima 
as shown in Fig. [7|d. By increasing ujq the system will 
probe thus a less shallow potential resulting in a stronger 
hardening of JIq than the harmonic one and leading to 
ctQ,o > 0.5. ii) for U 1.36, on the other hand, the shape 
of the lattice potential V{x) is more similar to the panel 
(f ) of Fig. Elwhere the zero point energy is lower than the 
unstable minima. Increasing loq the system would thus 
probe the flat potential regions around these two min- 
ima leading to an increase of the renormalized phonon 
frequency smaller than the harmonic one {aQg < 0.5). 
These effects become smaller as the Hubbard repulsion 
U increases since the minima shift higher in energy and 
the IE aQg approaches asymptotically 0.5 for U oo. 

A similar behavior occurs for A = 0.7 where however 
the system in the absence of electronic correlation is still 
metallic although close to the bipolaronic MIT. In this 
case however increasing U leads to a reduction of the IE 
on r^o until, for U 0.85, an^ becomes smaller than 0.5, 
although in a less marked way than in the A = 1.5 case 
because of the smaller electron-phonon coupling. The 
evidence of this anomalous IE aci^ < 0.5 even for small 
A suggests thus that this regime is not related to the 
closeness to the bipolaronic metal-insulator transition. It 
is worth to stress in addition that the evidence of these 
two regimes of correlation prompts out from the analysis 
of the IE acig while the renormalized phonon frequency 
itself r^o does not show any signature of this crossover. 



VI. CONCLUSIONS 

In this paper we have studied in details the isotope 
dependence of several physical quantities: the quasi- 
particle spectral weight, the renormalized phonon fre- 
quency, and the static component of local charge and spin 
susceptibilities in the metallic regime of the Hubbard- 
Holstein model within DMFT approximation. The work 
is carried out for the half-filled model in the metal- 
lic phase without broken symmetries showing anoma- 
lous behaviors approaching the pair and the Mott metal- 
insulator transitions. We have shown that the anomalies 
of isotope coefficients can be helpful in revealing differ- 
ent physical regimes. In particular, we have found that 
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in the adiabatic regime 7 = 0.1 we can identify two dif- 
ferent regimes according to the degree of correlation. For 
L/ < 1-5 the physical properties of the system are qualita- 
tively similar to those of the pure Holstein model where 
the approach to the bipolaronic metal-insulator transi- 
tion obtained within DMFT is reflected in a near diver- 
gence of the lEs characteristic of bipolaronic binding. In 
this regime we find a divergent az and a strong increase 
of close to the pair MIT while both and ap have 
a near divergence at Aph [ctvia +00 and dp —00). 
The second regime establishes for U > 1.5, where the 
bipolaronic MIT becomes of the first order. 

The change in the order of the transition can be char- 
acterized in an accurate way by means of our study of 
lEs. In the highly correlated regime, the isotope effects 
are indeed not diverging even when the first-order transi- 



tion is approached. The lack of any divergence as well as 
an analysis of the lattice potential suggest that the first- 
order character is more likely ascribed to an electronic 
mechanism than to a lattice-driven instability. In addi- 
tion we have shown that approaching the Mott-Hubbard 
MIT also gives rise to drastic anomalous lEs. 
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